An exactly solvable model of reversible adsorption on a disordered substrate 
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We consider the reversible adsorption of dimers on a regular lattice, where adsorption occurs on a 
finite fraction of sites selected randomly. By comparing this system to the pure system where all sites 
are available for adsorption, we show that when the activity goes to infinity, there exists a mapping 
between this model and the pure system at the same density. By examining the susceptibilities, we 
demonstrate that there is no mapping at finite activity. However, when the site density is small 
or moderate, this mapping exists up to second order in site density. We also propose and evaluate 
approximate approaches that may be applied to systems where no analytic result is known. 

PACS numbers: 



I. INTRODUCTION 

Over the years a significant and sustained research ef- 
fort has been directed at understanding the effect of dis- 
order in adsorption processes [H 01 ■ This is commensurate 
with the widespread occurrence of disorder in many dif- 
ferent kinds of adsorption, e.g. gases on solid surfaces and 
in porous media 3, biomolecules [J and colloidal par- 
ticles Q adsorbing from solution onto solid surfaces, 
and catalysis 0, ■ 

From a theoretical perspective, the challenge is how to 
incorporate disorder into statistical mechanical descrip- 
tions of the equilibrium and kinetic properties. In the 
absence of exact results, it is desirable to develop ap- 
proximate treatments, possibly by mapping the system 
containing disorder to a pure system with no disorder. 

One way to represent disorder is the random site 
surface (RSS) in which adsorption sites are uniformly 
and randomly distributed on a plane. The adsorbent 
molecules, represented by hard spheres, bind to these 
immobile sites. Steric exclusion is incorporated in the 
model in that a site is available for adsorption only if 
the nearest occupied site is at least one particle diameter 
away. The model is characterized by the dimensionless 
site density. 

This model was originally studied in the context of 
irreversible adsorption where it was shown that there ex- 
ists a mapping to an irreversible adsorption process on 
a continuous surface Q. The existence of this mapping, 
which is exact and valid in any dimension, means that 
from a knowledge of the amount adsorbed as a function 
of time in a continuous space we can calculate the amount 
adsorbed on the RSS surface of given site density at any 
time. Adamczyk and coworkers have successfully applied 
the model and its extensions (e.g. allowing for finite size 
adsorption sites) to the adsorption of latex spheres on 
mica surfaces 

More recently, Oleyar and Talbot studied the 
reversible version of the RSS model in which hard 
spheres adsorb and desorb from immobilized sites in a 
plane. Here the quantities of interest are the adsorption 
isotherm, i.e. the amount adsorbed as a function of the 



bulk phase activity of the solute and the structure of the 
adsorbed layer. Somewhat surprisingly, the theoretical 
description of the system in equilibrium (even when no 
phase transitions intervene) is more difficult than in the 
irreversible case since there appears to be no exact map- 
ping to the reversible adsorption on a continuous sur- 
face. Moreover, in one-dimension where the exact so- 
lution of the corresponding model without disorder, i.e. 
hard rods on a line is well known 11], attempts to solve 
the equilibrium RSS have so far proved unsuccessful (de- 
spite the fact that many one-dimensional statistical me- 
chanical models have exact solutions). 

It is the purpose of this article to study a simpler 
model, the adsorption of dimers on a one-dimensional 
lattice. Disorder is introduced by randomly eliminating 
a given fraction of the sites. The advantage of the lat- 
tice model is that essentially exact solutions are available 
for both the pure and disordered systems. This permits 
us to examine the existence of a possible mapping be- 
tween the two. We show that there is a mapping in the 
limit of small and large activities, but not in general. In 
the absence of a full mapping it is still useful to inves- 
tigate approximate approaches that permit an accurate 
description of the thermodynamics of the disordered sys- 
tem. We note that Oshanin et al. 0, Q studied a similar 
lattice model for the catalytic reaction A-\- A Q. 

We show that the introduction of an effective activity 
leads to highly accurate estimates of the thermodynamic 
properties of the disordered system. This methodology 
can be generalized to more realistic, off-lattice models. 



A. The pure model 

The system consists of adsorbed molecules which are 
in equilibrium with a bulk phase containing an adsor- 
bate at activity A. The molecules bind to sites of a one- 
dimensional lattice. Occupancy of one site excludes oc- 
cupancy of the nearest neighbor sites: See Fig. [T] We 
note that this model is isomorphic to the adsorption of 
dimers on the dual lattice [Q. For a system consisting 
of N adsorption sites the adsorbed phase can be formally 
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FIG. 1: Adsorption on a pure lattice. The filled circles 
resent the adsorption sites (A'^ = 30), while the open cii 
show the adsorbed molecules (n = 9). Steric exclusion 
vents the simultaneous occupancy of any two adjacent si 



FIG. 2: Adsorption on a diluted lattice. Out of a total of 30 
sites 14 (filled) are available for adsorption. These sites form 
clusters of lengths 2,2,4,1,1,3, and 1 



described with the grand canonical partition functioi 



N-l 
{zi=0,l} i=l 

here Zi is the site occupancy, A — exp(/3/i) is the activ- 
ity and free boundary conditions are imposed. We use 
the superscript * to denote the pure system, i.e. one 
in the absence of disorder. The exact calculation of the 
partition function is a simple exercise using the transfer 
matrix approach: See the appendix. 
The result can be expressed as 



5*(7V,A) = l+£; ( ^ )A" (2) 



where n^^ax = [{N + l)/2] is the maximum number of 
dimers that can be adsorbed on the lattice of size N {[x] 
represents the integer part of a;) 0, [1] . 

The average number of dimers adsorbed on a chain of 
N sites at an activity A is given by 



AAjv(A) = A 



d\iiE*{X,N) 
dX 



(3) 



and the fraction of occupied sites can be computed as 

J^n{X) 



p*iN,X) 



N 



which in the thermodynamic limit is given by 

2A 



P*(oo, A) 



(1 + V1 + 4A)V1 + 4A 



(4) 



(5) 



This has the expected behavior in the limits of small 
(Langmuir isotherm) and large (half of the sites are oc- 
cupied) activities. The susceptibility, or fluctuation in 
the number of adsorbed molecules, is given by 

<n^ > - <n>^ dp 
^^P^ N = ^dX 

after some calculation one obtains that, in the thermo- 
dynamic limit, 



X*(p) = p(l-p)(l-2p) 



(7) 



B. Model with disorder 

We now consider a diluted site model in which only 
a fraction of randomly selected sites are present. The 
grand partition function in this case is given by 

S(A, M) - E E n (1 - ^^^^+lV^V^+l)y"''' 

{n} z,=Q,i 1=1 

(8) 

where r]i — 0,1 denotes the absence or presence of an 
adsorption site i. The probability of finding an adsorbing 
site is given by 



P{ri) = ps<5,,,i + (1 - Ps)S,^fi. 



(9) 



There are no correlations between sites. A sample con- 
figuration is shown in Fig. O 

For a given site configuration {rji}, the number of clus- 
ters of exactly I contiguous adsorption sites is given by 



N 



ni 



{{"ni}) = ^ '7i-i)77i ■ • ■'y^+Ki - m+i+i)- (10) 



The mean number of clusters of size I is then given by 



ill) 



One can check that the sum rule characterizing the total 
number of occupied sites is verified, namely 



(12) 



1=1 



Note that the thermodynamic limit has been taken in 
the two equations: corrections for a finite system occurs 
for 1 — ps — 1/N and need not be included in these 
calculations. 

Since the adsorption site occupancies are quenched 
variables, the average over disorder is taken over the loga- 
rithm of the partition function, with the following result: 



ln(S(p„A)) = E^l^^*(^'0> 



(13) 



1=1 



where we have used the result for the partition function 
of a finite lattice of I connected sites (i.e. the pure model) 
with free boundary conditions. The mean density of ad- 
sorbed dimers is given by 



p{ps,X) = A 



51n(S(p„A)) 



dX 



E^M(A) (14) 



1=1 



3 



where 



M(A) = A 



dHE*il,X)) 
dX 



(15) 



is the average number of dimers on a (full) lattice of size 
/ with free boundary conditions. 

For small site density, one easily obtains 



A 



-Ps-2- 



l + A'^' (1 + A) (1 + 2 A) 
(3 + 2 A) A3 3 



Ps 



- 2 



(1 + 2A)(1 + 3A + A2)(1 + A)' 
A^ (3 A + 2) 



(1 + 4 A + 3 A2) (1 + 3 A + A2) (1 + 2 A) 



Ps^ + O (p/) 
(16) 



On the other hand, whatever the site density, when the 
activity A — > oo, A/i(A) ^ p for / = 2p and A// (A) p+l 
for Z = 2p + 1, so that one obtains 



piPs,Oo) = ( 



which gives 



p{Ps, oo) = 



1' 



(17) 



(18) 



a result identical to the case of hard rods onto a line 
with random sites[135- (Note, however, that ps cannot 
take any positive real value and is bounded by 1). 

In order to investigate the existence of a mapping be- 
tween this model and the pure model, we consider the 
thermodynamic quantities, in particular we focus on the 
susceptibility. 

In the presence of disorder, there are two kinds of sus- 
ceptibility: the connected susceptibility Xc{Ps, A) describ- 
ing the thermal fluctuations of the density of adsorbed 
particles and the disconnected susceptibility XdiPs, A) de- 
scribing the disorder-induced fluctuations of the density 
of adsorbed particles. 

The connected susceptibility is given by the thermo- 
dynamic relation (l3 | 



Xc(Ps, A) = A 



dpiPs,X) 
dX ■ 



(19) 



For large activities, the connected susceptibility goes to 
as 1/A. This result has a simple physical interpreta- 
tion. For a given realization of the disorder, the densest 
configuration is selected when A ^ oo. In this limit fluc- 
tuations are suppressed, and the connected susceptibility 
vanishes. 

The disconnected susceptibility is given by the second 
cumulant (in the average over the disorder) of the num- 
ber of adsorbed dimers for a given configuration of the 
disorder, Y.i^i({Vi}WiW- 

1 °° 

Xd{Ps,X) = ]^ H {nmF-mWWi{X)Afi'{^) (20) 



The different terms of the disorder average ninii can 
be obtained by sorting out configurations of overlapping 
and non-overlapping clusters: 

1. When / = I', there is a contribution when the 
two clusters are in the same location, and rTjjTiT = 
5u'N{l-p,fp[. 

2. When two clusters have a boundary in common, 
rrjnf — 2N{1 — ps)^p[^'' , the factor 2 coming from 
the two possibilities, right and left. 

3. When two clusters have overlapping sites, njn;/ = 
and the number of possibilities is equal to {I + I' + 

4. For all other configurations, there is no overlap be- 
tween clusters and nj/nj/ = nirTjJ'. this corresponds 
to N{N — I — I' — 3) configurations. 

After some calculation, one obtains that 
W>W -niW = N{1 - psfp[ 

[Sw + 2(1 - ps)p':+' -{l + l' + m- p.)Vi') • (21) 

In the large-activity limit, A ^ oo, one easily shows 
from Eqs. © and Q that 



OC 

E 

1=1 



Ps7V/(oo) = 



(l-p,)'(l + p.) 



(22) 



and 



From Eqs ^ a.nd ^ the following expres- 

sion of Xd{Ps, oo) now results: 



Xd(Ps,oo) = 



Pa(l - Ps) 
{^-Ps? ■ 



(24) 



Finally, by inverting Eg. pS)) . the disconnected suscepti- 
bility is obtained as 



XdiPs.co) = p{l - p){l - 2p) 



(25) 



which is identical to the susceptibility in a pure system 
at the same density (Eq.([7])). Since Xc(Ps, oo) = 0, this 
result establishes that the total fluctuations of the ad- 
sorbed number of dimers is the same in the disordered 
system at infinite activity and in the pure system at the 
same mean adsorbed density, p. This points to the ex- 
istence of an exact mapping between the two systems 
under these conditions. The mapping can actually be 
proven by extending the above considerations to the cor- 
relation functions. 

Is the mapping exact for finite activities? In order 
to answer this question, we compare x(p(Ps,A)) and 



4 



Xc(Ps,A) + Xd{Ps,^) for small values of ps- The first 
orders of the expansions read 



i2 



A 



Ps 



A + 1 
A^ (8A + 5) 2 



(1 + 2A) (A + 1)^ 
A^ (ISA^ + 69 A2 + 66A+ 17) 
(1 + 2A)(1 + 3A + A2)(A + 1)^ 
+ 0(p/)) 

where we have combined Eqs. [71 and [TBI and 

XciPsA) + XdifisA) = J^Ps 

_ A^ (8 A + 5) 2 
' (1 + 2 A) (A + 1)'^' 

(18A2 + 49A+ 17) 
(1 + 3A + A2)(A + 1)'(1 + 2A)^ 

+ 0{ps% 



Ps', 



(26) 



(27) 



It is noticeable that the expansions coincide at first 
and second order, but differ at third order and above. 
We speculate that this property is general in the sense 
that it is independent of the model (dimers, k-mers, hard 
rods,...) and of the dimension. In addition, one can see 
that the two expansions approach the same limit when 
the activity is very large. This means that they differ 
maximally for intermediate values of the activity. In any 
case, the above considerations prove that no mapping 
exists between the correlation functions of the disordered 
system at finite activity and those of the pure system at 
the same density, p{ps,X). 



II. APPROXIMATE SCHEMES 

The lattice model studied here is exceptional in that 
we know its exact solution, even when disorder is present. 
More often, we will have only a partial description, e.g. 
the cluster expansion that applies at low site density in 
the case of the RSS [l3|- In these cases, we seek approx- 
imate schemes that can provide an accurate description 
of the adsorption isotherms for the full range of site den- 
sities and activities. In this section we propose two such 
schemes and evaluate their accuracy by comparing with 
the exact results. 

One approach involves a partial resummation of the 
site density expansion of Eq.((T4]). Specifically, retaining 
explicitly the first term that appears at all orders we can 
write 



p(p.,A) = ^ - 



1=1 



-A 
1 + A 



+ Fiip,,X)\p[ (28) 



where Fi{ps, A) represents the remaining terms in the ex- 
act expansion and, consistent with these terms, has the 



0.4- 



0.3- 



0.2 



0.1 




FIG. 3: Density p as a function of activity A for ps = 
0.4,0.6,0.9 calculated via the exact formula, Ea. (|14[) (full 
curves), via Ea. (|30|l (dotted lines), and via Ea. (|36|l to the 
zeroth order (/(ps. A) = 1) (dashed curve) (see text). 



property that 



Fi{p,,X)^0 



(29) 



when A — > and A — > oo. 

The simplest approximation is to set Fi(X) = 0, result- 



ing in the approximate isotherm 



PiPs^X) 



Xps 



Xps + l + X 



(30) 



that, as expected, has the correct behavior in the limits 
of small and large activities. 

In order to highlight the deviations from the exact re- 
sults, we have plotted in Fig[3]the density p as a function 
of the activity A for various values of the site density ps ■ 
The approximate isotherms always overestimate the ad- 
sorbed density for all activities. The deviations increase 
with the site density and are most pronounced at inter- 
mediate activities. In order to understand the origin of 
this discrepancy, we have performed an asymptotic ex- 
pansion of Eq. (fT4)) at large activity. 



PiPsA) 



Ps 



1 p,(-p^ + 3p3 + 6) 

l + Ps 6A {l + ps)^il- Ps) 

and compared with that of Eq. |3D1 

Ps Ps 



PiPs,X) 



1 + Ps A(l+p,)2 



(31) 



(32) 



We note that when ps 0, Eq.dSIl) and Eq.(132]) coin- 
cide. This explains the very accurate description of the 
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FIG. 4: Effective activity Acff as a function of A for ps — 
0.9 calculated via the exact formula, Eq. (|14[) (full curve), via 
Ea. (|36|l to the zeroth order (dotted curve) and to the first- 
order, Ea. (|37|l (dotted-dashed curve). 



approximate isotherm in this limit. Conversely, when 
^ 1, one obtains 



A) 



1 + 6(1 -Ps)y 



(33) 



whereas no such combination in A(l — ps) appears in Eq. 
(1321). 

By including the leading term of the asymptotic be- 
havior, Eq. pTjl . the agreement with the exact result is 
only correct up to intermediate ps- It would be neces- 
sary to add additional terms when ps goes to 1 in order 
to have a good approximation of the isotherms. Such a 
procedure rapidly becomes very complicated and is not 
useful for more realistic models. The asymptotic analysis 
seems to indicate that the neglected terms of the series 
correspond to the existence of an asymptotic series whose 
variable is ((1 — ps)X)~^. 

An alternative approach uses an effective activity, Acff, 
instead of p{ps, A). The idea is to estimate the density in 
the disordered system with 



piPs,X) = P*(Aoff(A,p,.)) 



(34) 



where p* corresponds to the density at equilibrium given 
by Eq.([5]). Note that, since there is no exact mapping 
to a pure system for a finite activity A, even if we could 
find the function Acff (A, ps) that satisfies Eq. [MJ it would 
not give exact results for other thermodynamic quanti- 
ties. There is, however, merit in this approach because, 
as we have shown in previous sections, there is a map- 
ping in the limits of large and small activities. Thus, by 
using Eq. (PT|) with Eq. (O, one obtains the asymptotic 



expansion 

1 _ I 
Acff 



6 + 3ps~ Ps 
6psX 



o{\- 



(35) 



The coefficient of the leading term of this expansion is 
a function of p^ which behaves simply in the two limits 
Ps ^ and Ps — > 1- By combining with the low- activity 
expansion (1/Aeff = l/(Aps)), we propose the following 
interpolation scheme 



1 

Aeff 



(-1+Ps)' , /(A,ps) 



Ps 



Xps 



(36) 



where f{X,ps) is a [n,n\ Pade approximant. The 
isotherms are then calculated using Eq. (jM)) . As can 
be seen in Fig. [31 even to zeroth order (/(A,ps) = 1), 
this route gives a significantly better approximation than 
Eq. ipOl) . By matching the exact asymptotic behavior, 
Eq. ([33]) ■ as well as the second-order low-activity expan- 
sion, one obtains 



/(A,Ps 



1 + a(ps)A 
l + b{ps)X 



where 
a(ps 

and 



6 + 2 ps-^ + 9 ps^ - ps'^ 9 Ps - Ps 
6 - 3ps - Ps^ 



6 - 3ps - Ps"" 



(37) 



(38) 



(39) 



Since the zeroth order approximation is practically in- 
distinguishable from the exact result in the isotherm plot 
(Fig. [3]) , we highlight the differences between the two by 
plotting the effective activity versus the activity for high 
density ps in Fig. 3] The horizontal line corresponds to 
the asymptotic value, the full line to the exact result, the 
dotted line to the zeroth order, and the dashed line to 
the first-order. We note that this method converges very 
rapidly to the exact result, even at high ps, which is very 
encouraging for systems where no exact result is known. 

We can thus conclude from the present study that (i) 
There is a mapping between the adsorbed configurations 
of dimers on a diluted one-dimensional lattice and those 
on a full lattice at the same adsorbed density in the two 
limits of small and large activities; (ii) This mapping 
does not extend to finite activities. However, a success- 
ful approximation scheme is provided by considering a 
pseudo-mapping through the introduction of an effective 
activity. The scheme has been tested on the present one- 
dimensional lattice system for which we have also ob- 
tained exact expressions, and work is in progress to apply 
it to more realistic situations. 
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where Ai and A2 are the eigenvalues of the matrix T, 



FIG. 5: Adsorption isotherms on a pure chain of N sites 

{N = 2,4,8,16,32,64 top to bottom). The dashed Hne is 1 ± VTTlA 

the thermodynamic limit. Ai,2 = hr^- (A5) 



APPENDIX A: PURE MODEL 

To obtain an analytic expression for the partition 
function of the pure system we use the transfer matrix 
method: see, e.g. [ll|. Using this approach, Eq.dl]) can 
be rewritten as 

S*(A,iV) =Tr(T^yl) (Al) 

where 



For finite N the partition function is a polynomial in A, 
e.g. S*(A,1) = 1 + A, S*(A,2) = 1 + 2A, S*(A,3) - 
1 + 3A + A2, S*(A,4) = 1 + 4A + 3A2,... In the thermo- 
dynamic limit, only the largest eigenvalue contributes to 
the thermodynamic quantities and we have 

lim ln(S-(A,A^)),,/l + VTT4Ay 

N^oo N \ 2 J ^ ^ 

Some isotherms for different values of N are shown in 
Fig. [5l Note that the thermodynamic limit is approached 
rather slowly. 
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